Towards a realistic axion star 
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Abstract. In this work we estimate the radius and the mass of a self-gravitating system made of 
axions. The quantum axion field satisfies the Klein-Gordon equation in a curved space-time and 
the metric components of this space-time are solutions to the Einstein equations with a source 
term given by the vacuum expectation value of the energy-momentum operator constructed from 
the axion field. As a first step towards an axion star we consider the up to the 6 term in the 
axion potential expansion. We found that axion stars would have masses of the order of asteroids 
(~ 10 -10 M©) and radius of the order ~ few centimeters. 



1. Boson stars 

Boson stars (BS) are gravitationally bounded systems made of scalar particles. They 
were introduced for the very first time in 1968 by Kaup and later by Ruffini and 
Bonazzola [1]. The inclusion of a self-interacting term was done in [2] and the stability 
under general perturbations was studied in [3]. BS are fully characterized by the scalar 
field properties, i.e. the mass m of the scalar field and its potential V(0) = m 2 2 + 
A0 4 /2, where A is the self-interaction parameter. BS rise as solutions of the Einstein- 
Klein-Gordon equations 

G^v = 8ttG < T^ v > , (n-^L\<j)=0, (1) 

where □ = ( 1 / y/—g) dn [y/—gg^ v d v ] and the energy-momentum tensor is given by 

r MV = dp&dv* - \g»v (g aP d a <t>d p <t> + mW + |$ 4 ) . (2) 

Here < ... > denotes average over the ground state of the system of many particles. The 
field is second quantized. In this formalism, = + + 0~isan operator, where we 
have defined 0+ = ^+ mn R nl {r)Y l m {d , y^" 1 , <{,- = E/^J^(r)y£*(0, ¥ )e+' E "' and 
are the creation (annihilation) operators for a particle with angular momentum /, 
azimuthal momentum m and energy E n . From operator it is possible to construct the 
energy-momentum tensor (Q\Tn V \Q) by considering an state \Q) for which all the N 
particles are in the ground state (I = m = 0, n = 1). It is assumed an static, spherical 
symmetric space-time metric 

d s 2 = -B(r)dt 2 +A(r)dr 2 + r 2 d£l 2 . (3) 



By demanding regular solution at the origin and flatness at infinity, there exists a family 
of radial solutions characterized by the number of nodes in the radial coordinate. Those 



with zero nodes are called ground states. That solutions exhibit a maximum mass for a 
specific initial central field value C (O). It is known that boson stars with initial central 
field value 0(0) < C (O) are stable and those with 0(0) > C (O) are unstable. The magni- 
tude of the maximum mass of the boson star depends on the value of the self-interacting 
coupling. It was shown in [2] that maximum mass scales as M max = 0.22A 1//2 M 2 /m 

where M p is the Planck mass, A = XM 2 /4nm 2 is the adimensional self-interacting 
coupling and m the mass particle associated with the scalar field. 

2. The axion 

An axion is a scalar field originated at the Peccei-Quinn symmetry breaking phase 
transition [4] . It is one of the favorite candidates to account for the dark matter content 
of the universe. Its properties have been strongly constrained from astrophysical and 
cosmological considerations. The axion decay constant is allowed to be in the range 
10 10 GeV < f a < 10 12 GeV and its mass in the range 10~ 5 eV < m < 10~ 3 ev. At late 
stages of the evolution of the universe it acquires an non-vanishing potential energy 



density V(0) = m 2 f 2 
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3. Towards an axion star 

By doing and expansion of cos(0// a ) we obtain 
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The addition of the potential energy density (4) to the kinetic energy of axions give rise 
to a energy-momentum tensor of energy similar to eq. (2) except for a sign with the 
proper identification A = m 2 /6f 2 . It is natural to investigate if this energy density gives 
rise to standard boson stars as those studied in [2]. We proceed by following [1] with 
the potential energy eq. (4). The Einstein-Klein-Gordon equations (1) obtained with the 
metric (3) and (Q\T^ V \Q) computed as explained above are 
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where the prime denotes derivative with respect to x and we have done the definitions 

2 



x = rm, R\q 
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4. Results 



We have solved the system (5) for different values of A. As in the X \ 1 4 case, we found a 
maximum mass for specific c(0) = (7 C (0), but, the switch in the potential sign produces 
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FIGURE 1. Mass of a boson star as a function of the central value of the scalar field. Circles denote the 
critical point, that is, the maximum mass of a boson star. The stars with a central value to the right of this 
point are unstable and those to the left are stable. 



a significant change in the behavior of the maximum mass. The relation M max ~ A 1 / 2 is 
not satisfied anymore. Instead of increasing M max we found decreasing M ma x as can be 
easily observe in Fig. (1). We have extrapolated the maximum mass as a function of A 

»„ /n 0-844 1.396 0.0450 

M max = M max (A) = — + + . (6) 

Plugging the axion values of f a and m in A we found A ax i on ~ 10 13 — 10 17 and using this 
value in extrapolated M max eq. (6) we found 1O~ 12 M < M max (A ax j on ) < 10~ 10 M Q . The 
value <T C (0), where the maximum mass is reached, depends strongly on A. Again, by 
doing an extrapolation, we can estimate the critical density for and axion star. Another 
interesting result is that the radius doesn't depend on A but it depends on cr(0). Once 
we have estimated <7 C (0) for axion star we can then estimate the radius using the 
extrapolation. We found o c (A ax j on ) ~ 10~ 4 and the adimensional radius R(o c (A ax i on )) ~ 
10 2 , that is, R ~ few cm. for axion mass m ~ 10 -3 eV. 
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